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ABSTRACT 

The energy density of cosmic ray protons (CRp) in star-forming galaxies can 
be estimated from (i) 7r°-decay 7-ray emission, (ii) synchrotron radio emission, and 
(Hi) supernova rates. For most of the galaxies for which values of all these quanti- 
ties are known, the three methods yield consistent energy density estimates, ranging 
from O(10 _1 ) eVcm~ 3 in galaxies with low star-formation rates, to O(10 2 ) eVcm -3 in 
galaxies with high star-formation rates. The only cases for which the methods do not 
agree are the composite starburst/Seyfert 2 galaxy NGC1068, whose 7-ray emission 
originates in black-hole accretion rather than star formation, and the Small Magellanic 
Cloud, where the discrepancy between measured and estimated CRp energy density 
may be due to a small CR confinement volume. 
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1 INTRODUCTION 

Active star formation in galaxies leads to production of cos- 
mic ray protons and electrons (CRp, CRe) via the Fermi-I 
diffusive shock acceleration mechanism in supernova (SN) 
remnants (Fermi 1954; Ginzburg & Syrovatskii 1964; Bell 
1978; Protheroe & Clay 2004). 

Timescales of starburst (SB) activity in galaxies are 
comparable to galactic dynamical timescales, tsb ~ T dya ~ 
10 s yr. On the other hand, in a SB region the characteristic 
timescales for protons to gain energy (by the Fermi-I accel- 
eration process), and to lose it by collisions with thermally- 
distributed protons (leading to pion production and decay) 
and advection, are r+ ~ r_ ~ 10 5 yr. Indeed, assuming the 
Fermi-I process to be at work in a SN remnant, r+ = E/E = 
(AE/E^At = p- x At = (lO/^o.i) At ~ 10 5 yr, where 
At ~ 10 4 yr is a typical SN remnant lifetime, and (3o.i is the 
speed of the SN ejecta in units of 0.1c. The main energy- 
loss timescale is largely determined by the SN-driven outflow 
emanating from the SB region, i.e. t_ ~ r out ~ 10 5 yr; thus, 
r+ ~ r_ << tsb - 

The respective lengths of these timescales suggest that 
in galaxies, during typical episodes of star formation, a bal- 
ance can be achieved by CR between energy gains and losses. 
Under basic hydrostatic and virial equilibrium conditions in 
a galaxy, a minimum- energy configuration of magnetic field 
and CR may be attained; this is equivalent to having CR and 
magnetic fields in (approximate) energy equipartition (e.g., 



Longair 1994 - but see Beck & Krause 2005 for a critical 
view) . 

The equipartition assumption enables deduction of the 
CRp energy density in star forming galaxies, U p (the main 
contribution to the particle energy density), indirectly from 
the electron energy density (from radio synchrotron mea- 
surements) if a theoretically motivated injection proton-to- 
elecron (p/e) ratio is assumed. 

Another way to derive U p is based on measuring the 
GeV-TeV 7-ray emission, which is largely from CRp inter- 
actions with ambient gas protons, which produce neutral 
(•7r°) and charged pions; ty° decays into 7-rays. So 7-ray 
measurements provide a direct measurement of U p . Only 
recently have such measurements been possible for star- 
forming galaxies - and only for a handful of sources (see 
below) . 

Third, since energy-loss timescales are shorter than 
star-formation timescales, U p can be estimated from the ob- 
served rate of core-collapse SN and the deduced residency 
time of CRp in the insterstellar medium - once a fraction of 
SN kinetic energy that is channeled into particle acceleration 
has been assumed. 

This paper describes measures of galactic CR energy 
densities based on the methods outlined above. Although 
not strictly independent, these methods are based on very 
different observables to estimate U p : radio emission, HE 7- 
ray emission, and the rate of core-collapse SN. We find that 
the three methods yield consistent results on U p for a sample 
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of 10 galaxies with widely varying levels of star formation 
activity from quiescent to intense starbursts. These are the 
only galaxies for which 7-ray data (plus radio data and SN 
rates) are available (see Table 1). After reviewing the radio, 
7-ray, and SN methods (sect. 2,3,4), the corresponding val- 
ues of Up are presented in section 5. The results are discussed 
in section 6 and summarized in section 7. 



2 PARTICLES AND MAGNETIC FIELDS 

CRc populations consist of primary (directly accelerated) 
and secondary (produced via n ± decays) electrons. While 
the exact form of the steady-state CRe energy spectrum is 
not a single power law, at high energies the flattening of the 
spectrum due to Coulomb losses can be ignored, justifying 
the use of the approximate single power-law form. Let then 
the combined (primary plus secondary) CRe spectral density 
distribution be a single power-law 



iVe(7) = iV e ,o(l + X) 1' 



(1) 



where the electron Lorentz factor 7 is in the range 71 < 7 < 
72, N et o is a normalization factor of the primary electrons, \ 
is the secondary-to-primary electron ratio, and q > 2 is the 
spectral index (the equality holds in the strong-shock limit 
of the Fermi-I acceleration process). 

Ignoring the contribution of low-energy electrons with 
7 < 71, the electron energy density is U e = N Ci o (1 + 
X) m c c 2 7 1_9 d7, where 72 is an upper cutoff whose ex- 
act value is irrelevant in the the limit of interest 72 > > 71 . 
For q > 2 and 72 >> 71, 



N e ,o(l + X )m c c 2 1 2 1 - q /(q-2). 



(2) 



For a population of electrons described by Eq.(l), traversing 
a homogeneous magnetic field of strength B in a region with 
(a spherically equivalent) radius r s located at a distance d, 
and emitting a 1 GHz synchrotron radiation flux of /1 ghz 
Jy, standard theory yields 



AT e , (l + x) = 5.72 x 1CT 15 V a(qY 



250" 



(3) 



where a(q) is defined and tabulated in, e.g., Tucker (1975) 

and </> is defined as </> = ( oTifc)" 3 ^) 2 !^)- From 
Eqs.(2, 3) we then derive 



U. 



2 ' 96 xl0- 22 250V 



7i 



q+2 



(4) 

(1 + X) " («-2)o( 9 )" • 1 } 

In order to compute U e from Eq.(4) we need to specify 
71 and B. To do so we make the following assumptions: 
(i) We assume that the low-energy limit of the electron 
power-law spectrum, 71 , marks the transition (for decreasing 
energy) from Coulomb ((Gould 1972; Rephaeli 1979) to syn- 
chrotron losses. For an electron of energy 7, the synchrotron 
loss rate is 



"(-) 



dt ' s y n 



1.30 x 10~ 2 V 



B 



whereas the Coulomb loss rate is 



"(-) 
V At I 



dt ' coul 



1.2 x 10 _1 n e , th 



1.0 + 



In (7/n e , t h) 



75 



(5) 



.(6) 



(Rephaeli 1979; Sarazin 1999). We then simply assume that 



CRe lose energy via Coulomb scattering for 7 < 71 and via 
synchrotron cooling for 7 > 71. 

(ii) We assume equipartition between the energy densities 
of particles (CR electrons and protons) and magnetic fields, 
Up + U c — B 2 /8tt. This condition may actually be at- 
tained by strong coupling between all the relevant degrees 
of freedom in the SB region. In terms of the CR p/e energy 
density ratio, k, the particle-field equipartition condition is 
Up[\ + (1 + x)/«] = B 2 /8tt, so that 



7.44 x 10" 
1 + X 



7 1 2 -' ? 250 q/2 y> 
(9-2) a(q) 



(7) 



Inserting Eq.(7) into Eq.(5) we get (^) sy n oc . 
Once the value of n ejt h is specified (see Table 1), by equating 
Eqs.(5,6) we deduce 71. 

The secondary-to-primary electron ratio x, which ap- 
pears in Eqs.(7), depends on the injection p/e number ratio, 
r p / e = (mp/m c ) < - qi " i ~ 1 ^ 2 (see Appendix), and on the gas op- 
tical thickness to p-p interactions. Proton-proton (p-p) col- 
lisions lead to the the production of charged and neutral pi- 
ons. An electron is produced in the decay n~ — > /i~ + + v e 
followed by ^ — > e~ + 17, + V e . Given the branching ra- 
tios in p-p collisions, only a third of these collisions produce 
electrons. The mean free path of CR protons in a medium 
of density n p due to p-p interactions is A pp = (a pp np) _1 ; 
for protons with kinetic energy T ~fewTeV the cross sec- 
tion is a P p ~ 50 mb = 5 x 10 -26 cm 2 (Baltrusaitis et al. 
1984). For a typical SB ambient gas density n p ~ 150 cm" 3 , 
A pp ~ 43 kpc. The probability for a single CR proton to 
undergo a pp interaction in its 3D random walk through a 
region of radius r s ~ 0.25 kpc (also typical of SB nuclei) 
is then ^/3r s /X P p ~ 0.01. We then estimate that in typi- 
cal SB environments, characterized by relatively strong non- 
relativistic shocks (q^ = 2.2, see Appendix), the secondary 
to primary electron ratio is x — Xo a/3 (r s /A pp ) ~ 0.3. For 
quiescent environments, with typical values n p ~ 1 cm" 3 
and r s ~ 2.5 kpc, we find x — 0.03. The higher value found 
for starbursts is in approximate agreement with results of 
detailed numerical starburst models for energies <J lOMeV 
(plotted in, e.g., Paglione et al. 1996; Torres 2004; Domingo- 
Santamarfa & Torres 2005; De Cea et al. 2009; Rephaeli et 
al. 2010). 

To compute the CR p/e energy density, k (see Ap- 
pendix), we assume power-law spectra for the CR parti- 
cles: (i) the electron spectral index q c is deduced from the 
measured synchrotron radiation index a, according to the 
standard formula q c — 2a + 1; and (ii) the proton spec- 
tral index is assumed to be close to the injection value, 
q P ~ f/inj — 2.1 — 2.2, for the dense, CR-producing, SB 
environments hosted in the central regions of some galax- 
ies, and equal to the leaky-box value, q p = g; n j + S ~ 2.7 
(where S ~ 0.5 is the diffusion index) for more tenuous, CR- 
diffusing, quietly star forming galaxies. Values of k for our 
sample galaxies are reported in Section 5. 

Finally, we obtain an explicit expression for U p : 



8tt 



1 + 



1 + X 



7.44 x 10" 
1 + X 



r, K 1 

[ 1 + IT^ 



7 2 - 9 250 9/ ^ V 
X J (q-2) a{q) 
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Table 1. Star-forming galaxies: the data. 



Object 


4 11 


-PI 

' s 


f[3] 

■> 1 GHz 


a l4] 
NT 


[5] 

e, th 


r[6] 
TIR 


SFRW 


^ [8] 


M l9] 


r[10] 
L 1 


r [ll] 
' res 


Notes 




(Mpc) 


(kpc) 


(Jy) 




(cm -3 ) 


(erg/ s) 


(Mr.; /yr) 


(yr" 1 ) 




(crg/s) 

0/ 


(yr) 




Arp 220 


74.7 


0.25 


0.3 


0.65 


300 


45.75 


253 


3.5 


a 24+0-10 


< 42.25 


9.0E+3 


SB 


M82 


3.4 


0.26 


10.0 


0.71 


200 


44.26 


8.2 


0.25 


Q o 7 +0.09 


40 21+ - 10 
^ u - zi -0.13 


2.6E+3 


SB 


NGC 253 


2.5 


0.20 


5.6 


0.75 


400 


44.23 


7.7 


0.12 


9 20+ - 10 


39.76±°-i| 


2.0E+4 


SB 


Milky Way 




4.4 






0.01 


43.75 


2.5 


0.02 




38.91±°;^ 


2.7E+7 


quiescent 


M31 


0.78 


5.17 


4.8 


0.88 


0.01 


42.98 


0.43 


0.01 


9.881°;" 


38.661°:°? 


4.0E+7 


quiescent 


M33 


0.85 


2.79 


3.30 


0.95 


0.03 


42.68 


0.22 


0.003 


q oc+0.13 


< 38.54 


2.6E+7 


quiescent 


LMC 


0.049 


3.0 


285.0 


0.84 


0.01 


42.45 


0.16 


0.002 


8.86t°;» 


37 67+ - 05 
d '-°'-0.05 


1.0E+7 


quiescent 


SMC 


0.061 


1.53 


45.3 


0.85 


0.01 


41.45 


0.01 


0.001 


o fifi +0.03 
8 - bb -0.06 


37 nA+O-ll 


4.0E+7 


quiescent 


NGC 4945 


3.7 


0.22 


5.5 


0.57 


300 


44.02 


4.7 


0.1-0.5 


9 64+ - 10 


4O.30t°;£ 


4.6E+4 


SB+Sy2 


NGC 1068 


16.7 


1.18 


6.6 


0.75 


300 


45.05 


50 


0.2-0.4 


q 71+O.H 
y - ' -"--0.19 


41 32+ 015 


1.0E+6 


SB+Sy2 



M Distance (from Abdo et al. 2011). 

PI Effective radius of star-forming region. See text. Data are from Persic & Rephaeli 2010 and refs. therein (Arp 220, M82, NGC 253), 
Beck & Grave 1982 (M31), Tabatabaei ct al. 2007 (M33), Weinberg & Nikolaev 2001 (LMC), Wilke et al. 2003 (SMC), Moorwood & 
Oliva 1994 (NGC 4945), Spinoglio et al. 2005 (NGC 1068). 

Ml GHz flux density. Data arc from Persic & Rephaeli 2010 and refs. therein (Arp 220, M82, NGC 253)), Beck & Grave 1982 (M31), 
Tabatabaei et al. 2007 (M33), Klein ct al. 1989 (LMC), Hayncs ct al. 1991 (SMC), Elmouttie ct al. 1997 (NGC 4945), Ruhr ct al. 1981 
(NGC 1068). 

' 4 1 Non-thermal spectral radio index. Data are from Persic & Rephaeli 2010 and refs. therein (Arp 220, M82, NGC 253), Beck & Grave 
1982 (M31), Tabatabaei et al. 2007 (M33), Klein et al. 1989 (LMC), Haynes et al. 1991 (SMC), Elmouttie et al. 1997 (NGC 4945), 
Ruhr et al. 1981 (NGC 1068). 

I 5 1 Thermal electron density. Data arc from Roy ct al. 2010 (Arp 220), Petuchowski ct al. 1994 (M82), Kewley et al. 2000 and Carral et 
al. 1994 (NGC 253), Cox 2005 (Milky Way), Beck 2000 (M31), Jabatabaei et al. 2008 (M33), Points et al. 2001 (LMC), Sasaki ct al. 
2002 (SMC), Spoon ct al. 2000 (NGC 4945), Kewley et al. 2000 (NGC 1068). 
^Total IR [i.e., (8 - 1000)/im] luminosity, in log (from Abdo ct al. 2011). 
PlStar formation rate, from SFR = L m /(2.2 X 10 43 crg/s) (Kcnnicutt 1998). 

M Core-collapse SN rate. Data are from Persic & Rephaeli 2010 and references therein (Arp 220, M82, NGC 253), Dichl ct al. 2006 
(Milky Way), van den Bergh & Tammann 1991 (M31, M33, SMC, LMC; see also Pavlidou & Fields 2001), Lcnain ct al. 2010 and 
references therein (NGC 4945, NGC 1068). For NGC 1068 we also computed an upper limit to the SN rate (iajn < 0.39) using Mannucci 
et al.'s (2003) formula u SN = (2.4 ±0.1) X 10- 2 [L PIR /(10 10 L Q )] yr" 1 , being / FIR = 1.26 X 10" 11 (2.58 f 60 + f 100 ) ergem^s- 1 (see 
Helou et al. 1988) with / 60 ~ 190 Jy and /100 ^ 277 Jy. 

I 9 ' Gas mass (neutral plus molecular hydrogen: Mjji + Mjj 2 ), in log- Data are from: Torres 2004 for Arp 220; Abdo et al. 2010a for M82, 
NGC 253, and the Milky Way; Abdo et al. 2010b for M31 and M33; Abdo ct al. 2010c for the LMC; Abdo ct al. 2010d for the SMC; 
and Lcnain ct al. 2010 for NGC 4945 and NGC 1068. 

[10] High-energy (>100MeV) 7-ray luminosity, in log (from Abdo et al. 2011). 
[ n lCRp residency time. 



Using Eq.(8), numerical values of U p can be obtained from 
the relevant observational quantities for our sample galaxies 
(see Table 1); these values are reported in Table 2. 



3 ESTIMATING U P FROM 7-RAY EMISSION 

In this section we will review some basic features of the mod- 
eling of 7-ray emission from star-forming galaxies, and the 
status of HE/VHE 7-ray observations. Detections of several 
such galaxies have enabled measurements of U p values either 
in SB cores or throughout the disk. 

In most SB galaxies, such as the two nearby ones M 82 
and NGC 253, the central SB region (also called the source 
region) with a radius of ~ 200 — 300 pc, is identified as the 
main site of particle acceleration. Here, the injection parti- 



cle spectrum is assumed to have the non-relativistic strong- 
shock index q = 2. A theoretical N p /N c ratio, predicted from 
charge neutrality of the injected CR, is likely to hold in this 
source region - as is also the assumption of cquipartition. A 
measured radio index a ~ 0.75 in the source region implies 
q = 2 a + 1 ~ 2.5 there. This indicates a substantial steep- 
ening of the CRe spectrum from the injection value qo = 2 
due to diffusion (D oc 7~ 4 ) effects that cause the steady- 
state particle spectral index to be qo + S above some break 
energy. As an example, the steady-state electron and proton 
spectra in the SB region of NGC 253 are shown in Fig. 1. At 
low energies both spectra are flatter, whereas at E >> 1 
GeV the stronger electron losses result in steeper electron 
spectra. 

Adopting the convection-diffusion model for energetic 
electron and proton propagation and accounting for all the 
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emission from p-p-induced n° decay is 



r M _ 

L r J >E — 



Figure 1. Properties of the emitting particles in the central SB 
region of NGC 253 (Rephaeli et al. 2010): steady-state spectra 
of primary (solid line) and secondary (dot-dashed line) electrons, 
and of protons (dashed line). 




Figure 2. Properties of the emitted radiation in the central SB 
region of NGC 253 (Rephaeli et al. 2010): radiative yields from 
electron Compton scattering off the FIR radiation field (dotted 
line), electron bremsstrahlung off ambient protons (dashed line), 
7T° decay following p-p collisions (dashed-dotted line), and their 
sum (solid line). 



relevant hadronic and leptonic processes, the steady-state 
energy distributions of these particles, in both the (active) 
SB nucleus and the (passive) disk of these galaxies, can be 
determined with a detailed numerical treatment once the gas 
distribution is known (e.g., Torres 2004; Persic, Rephaeli & 
Arieli 2008; Rephaeli, Arieli & Persic 2010). 

The relevant electron energy loss processes are 
bremsstrahlung, Compton, and synchrotron, whereas for 
protons the main losses are 7-ray emission from ty° decay 
following p-p collisions (see Fig. 2); bremsstrahlung losses 
dominate at lower energies, whereas 7r°-decay losses domi- 
nate at higher energies. In the GeV-TeV region, emission is 
mainly from p-p-induced 7r° decay. 

The procedure is similar when star formation is not un- 
dergoing a burst confined to the nuclear region but occurs 
throughout the whole disk. 

For a source with gas number density n gas , proton en- 
ergy density U p , and volume V, the integrated hadronic 



fl>L n gas Up dV s 1 



(9) 



with the integral emissivity g!?' in units of photons 1 (H- 
atom)~ 1 (eV/cm 3 )" 1 (Drury et al. 1994). Therefore U p can 
be determined, once L> e and n gas (r) are observationally 
known and the particle steady-state energy distributions 
have been numerically worked out by solving the convection- 
diffusion model for CRe and CRp propagation. By its very 
nature, this is a direct measurement of Up. 

The two local SB galaxies M 82 and NGC 253 are the 
only non-AGN extragalactic sources that have so far been 
detected in both the GeV (Abdo et al. 2010a) and TeV 
(Acciari et al. 2009; Acero et al. 2009) regions. The mea- 
sured fluxes and spectra of both galaxies in the two bands 
agree with predictions of recent numerical models in which 
U p = O(10 2 )eVcm" 3 in the SB nucleus. The highest-SFR 
galaxy in the nearby universe, Arp220, was undetected by 
MAGIC (Albert et al. 2007). 

HE 7-ray detections were obtained with the Large Area 
Telescope (LAT) on board the orbiting Fermi telescope for 
a number of low-SFR galaxies: (i) the Andromeda galaxy 
M31 (Abdo et al. 2010b), with Up ~ 0.35 eV cm" 3 ; (ii) 
the Large Magellanic Cloud (LMC) whose average spec- 
trum, either with or without the bright star-forming region 
of 30Doradus, yields Up ~ 0.2 - 0.3eVcm -3 (Abdo et al. 
2010c); (Hi) the Small Magellanic Cloud (SMC), whose in- 
tegrated spectrum indicates that Up ~ 0.15 eVcm~ 3 (Abdo 
et al. 2010d); and (iv) the composite SB/Seyfert-2 galaxies 
NGC 1068 and NGC 4945 (Lenain et al. 2010). Only flux up- 
per limits exist for the Triangulum galaxy M 33 (Abdo et al. 
2010b). 

For the Milky Way, the modeling of the Galactic dif- 
fuse HE emission along the lines outlined above suggests an 
average U p ~ leVcm -3 (Strong et al. 2010; Ackermann et 
al. 2011). This compares with t/ p ~leVcm -3 measured at 
Earth (e.g., Webber 1987), and with (6±3)eVcm -3 in the 
~200 pc region of the Galactic Center, as inferred from the 
measured VHE 7-ray emission (based on HESS data: Aha- 
ronian et al. 2006). 

The Up values estimated from analysis of GeV- TeV ob- 
servations are reported in Table 2. 



4 ESTIMATING Up FROM SN RATES 

Combining the SN frequency with the residency timescale, 
Tr CS , of CR protons that give rise to TeV emission in the 
star-forming region, and assuming a bona-fide value of the 
energy that goes into accelerating CR particles per SN event, 
we can obtain a second estimate of Up that stems from the 
link between CR and core-collapse SN. 

The value of r ros is determined from those of two other 
timescales: 

(i) energy-loss timescale for p-p interactions, r pp = 
(o"ppCn p ) -1 that, for protons with kinetic energy 10 TeV 



for which a v 



50 mb, can be written as 



2 x 10 5 {- 

v - 



100 cm" 3 ' 
(ii) CRp advection timescale, r, 



yr; (10) 
out, that characterizes the 
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removal of CR out of the disk mid-plane region in a fast 
(vout ~ 2500 km s" 1 for M82: Chevalier & Clegg 1985; 
Seaquist & Odegard 1991; Strickland & Heckman 2009) SB- 
driven wind, which for a homogeneous distribution of SNe 
within the SB nucleus of radius r s is 

r out = 3 x 10 4 (-^—) ( ^^—y 1 yr . (11) 

V 0.3kpc' V 2500kms- 1 '' y V ' 

(Both timescales are normalized to typical SB conditions.) 
Thus, 

T"res = fp'pW +T ut(r s ,iw)- (12) 

During r rcs , a number fsNfres of SN explode and de- 
posit the kinetic energy of their ejecta, E c j = 10 51 erg per 
SN (Woosley & Weaver 1995), into the interstellar medium. 
Arguments based on the CR energy budget in the Galaxy 
and SN statistics suggest that a fraction rj ~ 0.05 — 0.1 of 
this energy is available for accelerating particles (e.g., Hig- 
don et al. 1998; Tatischeff 2008). Accordingly, we express 
the CRp energy density as 

p 0.3yr-! 3xl0 4 yr 0.05 10 51 erg ^O.SkpC 1 

eV cm~ 3 . (13) 

5 VALUES OF U P IN GALAXIES 

In SB galaxies star formation activity is intense in a rel- 
atively small nuclear region, in addition to quiescent, low- 
intensity star formation throughout the rest of the galactic 
disk (which is characteristic of non-SB spiral galaxies). This 
is the case, for example, in the nearby M 82 and NGC 253. 
The star formation rate in the nuclear SB region can be suf- 
ficiently high that the main contribution to the CR density 
comes from acceleration (by SN shocks) in this region, with 
radial extent r s . The rest of the disk can then be viewed as 
target to the diffusing CR. In non-SB spiral galaxies, e.g. the 
Milky Way and the Triangulum galaxy M 33, star formation 
proceeds throughout the disk, where both CR sources and 
target gas clouds are distributed; in this case, r s is essentially 
the disk radius. 

For the sample of galaxies detected in 7-rays with 
Fermi/L AT (see Table 1), values of U p were deduced directly 
by modelling the observed 7-ray emission. Using Eq.(8) and 
Eq.(13), we now provide estimates of U p from measured ra- 
dio emission from CR electrons and the statistics of SN ex- 
plosions. 

Arp 220. The starburst activity takes place in a central 
molecular-gas disk of radius 480 pc and thickness 90 pc that 
accounts for SFR~120 Mo yr - , and embeds (at a galac- 
tocentric distance of 200 pc) two extreme starburst nuclei 
of radii 68 pc and 110 pc that account for SFR= 50 and 
35M yr _1 , respectively. This scheme, derived from Torres 
(2004) and using Eq. 4 of Kennicutt (1998), is consistent 
with the galaxy-wide IR-inferred SFR~ 225M yr _1 . For 
the disk alone we assume r s ~ 0.25 kpc, /1 ghz — 0.3 Jy, and 
a ~ 0.65, which implies q c = 2.30, and hence a ~ 0.09; 
assuming q p = 2.2 leads to k ~ 38. From Eqs.(6) and 
(5), and the value of n e ,th reported in Table 1, we derive 
71 = 2650. From Eq.(8) we then obtain B « 207 /iG and 
Up « 1027 eVcm~ 3 (radio method). In Eq.(ll) we use 
0.09 kpc, the thickness of the starburst disk perpendicular to 



which the wind breaks out, to compute r a( jv — 9 x 10 yr (be- 
ing n p ~ 10 2 cm -3 ; and using the latter and ^sn = 3.5 yr -1 
in Eq.(13), we get U p « 515eVcm -3 (SN method). 

M82. Following Persic et al. (2008; and references therein) 
we take the central SB to be a region with a radius of 300 pc 
and height of 200 pc, hence r B — 260 pc, whose non-thermal 
radio emission has /ighz = 10 Jy and a ~ 0.71. The latter 
implies q ~ 2.42, that yields a ~ 0.09; assuming q p — 2.2, we 
then compute k ~ 60. From Eqs.(6) and (5), and the value 
of n e , th reported in Table 1, we derive 71 = 4724. From 
Eq.(8) we then obtain B « 91 fiG and U p « 201eVcm~ 3 
in the SB nucleus (radio method). Because ^sn = 0.25 yr -1 
and T rcs ~ Tout — 3 x 10 4 yr (taking n p ~ 10 2 cm -3 in the 
SB nucleus, and v ou t ~ 2500 kms~l; see Persic & Rephaeli 
2010), from Eq.(13) we get U p « 95eVcm^ 3 (SN method); 

NGC 253. Following Rephaeli et al. (2010, and references 
therein) we assume r s = 200 pc, /ighz = 5.6 Jy, and a ~ 
0.75. The latter implies q ~ 2.50, hence a — 0.0852; assum- 
ing q p — 2.2, we get k ~ 76. From Eqs.(6) and (5), and the 
value of n ej th reported in Table 1, we derive 71 = 7886. From 
Eq.(8) we then obtain B « 77 /iG and U p « 145eVcm~ 3 in 
the starburst nucleus (radio method). Being issn = 0.12 yr -1 
and r rcs ~ Tout ~ 3 x 10 4 yr (being n p ~ 10 2 cm -3 , in the 
starburst nucleus, and v out ~ 2500 kms~l; see Persic & 
Rephaeli 2010), we get U p w 75eVcm~ 3 (SN method); 

Milky Way. Measurements of Galactic CR indicate U p ~ 
1 eV at the Sun's position. The Galactic CRp flux mea- 
sured locally is the result of the superposition of parti- 
cles, accelerated in several sites scattered throughout the 
Galaxy, that have diffused out into the disk. Accordingly, 
we view the entire Galaxy as the region where SN oc- 
cur, and adopt 15 kpc and a thickness of 0.5 kpc, so that 
r s = 4.4 kpc, r rcs ~ 2 x 10 7 yr (with n p ~ 1cm -3 ), and 
IA3N = 0.02 yr~\ leading to U p w leVcm -3 (SN method). 
Alternatively, if we directly compare our SN-based estimate 
to available measurements of the innermost Galactic re- 
gion, we need to estimate the relevant quantities in a nu- 
clear region with r s = 0.2 kpc. The most prominent large- 
scale distribution of HI gas in the Galaxy is an exponential 
disk with R d — 3.75 kpc coplanar to the stellar disk with 
rim = 1 cm" 3 in the solar neighborhood, i.e. at a Galacto- 
centric distance Ro = 8 kpc (Dehnen & Binney 1998): this 
implies that the average density within r B is rim — 8 cm -3 , 
hence r pp = 2.5 x 10 6 yr there. The full disk Galactic SN 
rate, i^sn = 0.02 yr -1 (Diehl et al. 2006), should similarly 
be rescaled to the r s region: In an exponential disk with 
i?d = 2.5 kpc (Dehnen & Binney 1998), r B contains 0.3% of 
the total mass; consequently, usk = 0.6 x 10 _4 yr _1 , leading 
to U p w 5eVcm~ 3 (SN method). 

M 31. With an average value of the spectral index, a = 0.88, 
and the (nonthermal) flux (at 2.7 GHz) reported by Beck 
& Grave (1982), we deduce /ighz = 4.8 Jy, averaged over 
19.2 kpc radius. The latter corresponds to ~3.34 exponential 
length scales (Rd = 5.75 kpc, averaging VRI data from Hou 
et al. 2009), i.e. the corresponding volume encloses ~85% of 
the mass of the exponential disk, hence most of the stellar 
distribution. With an assumed thickness of 0.5 kpc, we get 
r s — 5.17 kpc. The measured spectral radio index implies q = 
2.76, and a ~ 0.08; assuming q p — 2.7, we get k ~ 8. From 
Eqs.(6) and (5), and the value of n Si th reported in Table 1, 
we derive 71 = 1169, and from Eq.(8) we then obtain B w 
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2.6 /iG and Up ~ 0.15 eVcm (radio method). Taking the 
whole disk as a site for SN explosions, the implied average 
gas density is n p ~ 0.5 cm~ 3 , hence t icb ~ r pp ~ 3.7 x 10 7 yr. 
Eq.(13) then yields Up « 0.2eVcm~ 3 (SN method). 

M33. Star formation - traced by OB associations and HII 
regions (e.g., Bastian et al. 2007) - and SN occur through- 
out the disk. From Tabatabaei et al. (2007) we derive a non- 
thermal radio flux /ighz = 3.3 Jy and a — 0.95, averaged 
over 7.6 kpc radius. The latter corresponds to 5.3 exponen- 
tial lengthscales (R d ~ 1.43 kpc, from Regan & Vogel 1994), 
enclosing 97% of the mass of an exponential disk, hence vir- 
tually the whole stellar distribution. For an assumed thick- 
ness of 0.5 kpc, we get r s = 2.79 kpc. The measured spectral 
index implies q = 2.90, and a ~ 0.08; assuming q p — 2.7, 
we get k ~ 11.2. From Eqs.(6) and (5), and the value of 
n e ,th reported in Table 1, we derive 71 = 1287. From Eq.(8) 
we then obtain B « 4.1 /iG and U p ~ 0.38 eV cm -3 (radio 
method). Given the gas mass reported in Table 1, the aver- 



Table 2. Star- forming galaxies: CRp energy dcnsitics+. 



age gas density is n p ~ 1 cm , hence T r( 



2x 10'yr. 



Inserting the relevant quantities (see Table 1) into Eq.(13), 
we then derive U p ~ 0.7eVcm~ 3 (SN method). 
LMC. This satellite of the Milky Way can be modeled 
as a truncated disk/spheroid with r t ~ 3 kpc whose half- 
thickness is also ~3kpc (Weinberg &: Nikolaev 2001), so we 
use r s = 3 kpc in Eqs.(13,8). The measured spectral radio 
index implies q — 2.68, and a ~ 0.08; assuming q p = 2.7, we 
get k ~ 6.7. From Eqs.(6) and (5), and the value of n e , th 
reported in Table 1, we derive 71 = 1189. Eq.(8) then yields 
B « 2.6 fiG and U p « 0.14eVcm~ 3 (radio method). There 
is no mass outflow in the LMC, hence r ros ~ r pp w 10 7 yr 
(using an average gas density of n p « 2 cm -3 .) Eq.(13) then 
yields U p « 0.2eVcm~ 3 (SN method). 

SMC. This other Milky Way satellite can be modeled as 
a bar with an area ~ 2.5 x 1.5 kpc (Wilke et al. 2003) and 
l.o. s. extent of 4 kpc (following Abdo et al. 2010d), so we 
use r s = 1.53 kpc in Eqs.(13,8). The measured radio spectral 
index implies q — 2.70, and a ~ 0.08; assuming q p — 2.7, 
we get k ~ 20. From Eqs.(6) and (5), and the value of n e , th 
reported in Table 1, we derive 71 = 754. From Eq.(8) we 
obtain B w 4/iG and U p w 0.39 eV cm -3 (radio method). 
The SMC has a (galaxy-wide) SN rate of v S n = 10~ 3 yr~\ 
There is no mass outflow from the SMC, so r rcs ~ r pp ~ 
1.4 x 10 7 yr- (It is n p ~ 1.4cm .) From Eq.(13) we then 
obtain U p « leVcm~ 3 (SN method). 

NGC 4945. This almost edge-on galaxy hosts a Seyfert- 
2 nucleus (deduced from its X-ray variability, Iwasawa et 
al. 1993). However, essentially all its IR radiation arises 
from a central 33" x 19" (r s = 0.22 kpc) molecular complex 
whose IR colors are typical of SB activity (see Moorwood 
& Oliva 1994). The measured nuclear radio emission im- 
plies q = 2.14, and a ~ 0.10; assuming q p = 2.2, we get 
k ~ 15.7. From Eqs.(6) and (5), and the value of n e , th 
reported in Table 1, we derive 71 = 5689. Assuming the 
measured radio flux to be only related to the SB, from 
Eq.(8) we then get B « 94 fxG and U p « 201 eV cm" 3 (radio 
method). It follows from psn ~ 0.2 yr -1 , and the existence 
of a galactic wind with « out ~ 1200 kms -1 (Chen & Huang 
1997) that r rcs ~ r out ~ 4.6 x 10 4 yr; Eq.(13) then yields 
Up « 220eVcm" 3 (SN method); 

NGC 1068. This galaxy hosts a prototypical Seyfert-2 nu- 
cleus (e.g., Wilson & Ulvestad 1982). The circumnuclear SB 
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(a) Acciari et al. 2009 (see also Persic et al. 2008 and De Cea 
et al. 2009). (b) Accro et al. 2009 (see also Paglione et al. 1996, 
Domingo-Santaman'a & Torres 2005, and Rephaeli et al. 2010). 
(c) Abdo et al. 2010a. (d) Strong et al. 2010. (e) Aharonian et al. 
2006. (f) Abdo et al. 2010b, with Drury et al. 1994 in the case of 
M33. (g) Abdo et al. 2010c. (h) Abdo et al. 2010d. (i) Lenain et 
al. 2010. (j) Webber 1987. 

region can be modeled as a spherical shell with external ra- 
dius of 1.5 kpc and thickness 0.3 kpc, and mass 3.4 x 1O 9 M0 
(Spinoglio et al. 2005). The effective radius is then r s — 
1.2 kpc, and the average density is n p ~ 20cm~ 3 . The mea- 
sured spectral index implies q ~ 2.50, hence a = 0.0852 
- assuming q p — 2.2, we get k ~ 76. From Eqs.(6) and 
(5), and the value of n e ,th reported in Table 1, we derive 
71 = 10161, and from Eq.(8) we then compute B w 52 piG 
and U p w 65eVcm~ 3 (radio method). The CRp residency 
time is r ros ~ r pp ~ 10 6 yr. (The outgoing wind observed in 
NGC 1068 does not emanate from the SB region, but from 
close to the central AGN, see Krolik & Begelman 1986.) As- 
suming all the observed far-IR emission originates from the 
SB region, yields usk = 0.39 yr -1 there (see Table 1). Insert- 
ing the relevant quantities into Eq.(13) we get B w 52 /j,G 



and U p w 65eVcm- 3 (SN method). 



5.1 Uncertainties 

The values of U p listed in Table 2 do not include substantial 
observational and modeling uncertainties. In this section we 
attempt to estimate the level of precision with which Up was 
determined based on only limited information on the errors 
in the various observational parameters. 



5.1.1 Radio Method 

The quantities in Eq.(8) are usually well determined for our 
sample galaxies, except for the p/e energy density ratio k, for 
which a CRp spectral index, q p , must be assumed. Given its 
possible values discussed in the Appendix (i.e., q p ~ 2.1 — 2.2 
in SB regions, and q p ~ 2.1 — 2.2 for quiescent galaxies), the 
uncertainty on the CRp spectral index, 5q p ~ 0.1, translates 
to a factor of ~2 uncertainty on k [see Eq.(24)], i.e. typically 
an uncertainty of ~50% on Up as deduced from Eq.(8). 
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We emphasize that the use of k here is slightly different 
from that of Persic & Rephaeli (2010), where we assumed 
q P = q c — 2a + 1 in all galaxies, whereas here we appro- 
priately drop such equality and assume a physically moti- 
vated index, separately for SB (q p gi n j) and quiescent 
(q p ginj + S) galaxies. Numerically the values of q p as- 
sumed here are not that different from the measured q c . 

Another improvement in the present radio-based deriva- 
tion of U p relates to the electron low-energy cutoff, 71. In 
Persic & Rephaeli (2010) we assumed 71 = 10 3 , whereas here 
we calculate 71 self-consistently by setting it equal to the 
energy at which synchrotron and Coulomb losses are equal. 
Effectively then, 71 is the low-energy limit of the electron 
PL spectrum (see Eq.l). 

5.1.2 Supernova Method 

A precise measurement of the actual rate of core-collapse 
SNe is, of course, crucial in these estimates. This is quite 
difficult given the typically heavy optical extinction (e.g., in 
SB nuclei). For example, from Table 2 notice that for our SB 
nuclei the SN-based U p are systematically lower, by a factor 
of 2, than the corresponding radio and 7 estimate, whereas 
there is no such discrepancy in the quiescent sample. This 
mismatch does suggest the possibility of missing about half 
of the SN remnants hosted in dense, dusty environments. In 
addition, SN rates, even when quoted for individual galax- 
ies, are often only statistical in nature and based on the 
morphology and stellar mass of the respective galaxies (e.g., 
Mannucci et al. 2005), especially so at low luminosities where 
the SN rate is very low. Finally, radio counts of SN remnants 
require information on their ages in order to be turned into 
actual SN rates. For our sample galaxies, published observa- 
tional results suggest that z^sn are known to within a factor 
~2. 

Also from Eq.(13) we see that uncertainties on U p can 
also be due to uncertainties in r ros . The latter mostly arise 
from the uncertainty in the fast wind velocity, which in turn 
is probably known to within a factor of ~2. Estimates of 
the SN energy that is channeled to CR, ~ (0.5 - 1) 10 50 erg, 
i.e. ~5% — 10% of the total kinetic energy of the SN ejecta 
(Woosley & Weaver 1995), agree within a factor ,$ 2 (Lin- 
genfelter et al. 1998; Higdon et al. 1998). (Here we assume 
that all core-collapse SN share a universal CR acceleration 
efficiency.) 

Finally, for all our sample galaxies, the effective star- 
forming radii r B are deduced from high-resolution optical 
and radio data, so their values are relatively precise, with 
typical uncertainties of a few 10%. 

We conclude that all the galaxy quantities in Table 1 
relevant to Eq.(13) are quoted in the literature as observa- 
tionally precise to within a factor of ,$ 2, so that U p can 
only be estimated from Eq.(13) to within a factor ~4. 



6 DISCUSSION 

The three methods discussed here are not independent. 
The 7-ray method and the radio method are coupled 
through the p/e ratio at injection, through the secondary- 
to-primary electron ratio, and through the imposed condi- 
tion of particle-field equipartition. The SN method is not 



independent of the 7-ray method either, because both de- 
pend on the residency time of CRp in the emission region 
- although, unlike the 7-ray and radio methods, it does not 
depend on the particle radiative yield but on the statistics 
of core-collapse SN. Also, the 7-ray, radio, and SN methods 
are not on equal footing. By these methods the value of U p 
is either measured, inferred, or estimated, respectively. This, 
because: (i) 7r°-decay 7-ray emission is the most robust mea- 
sure of U p when the distribution of target gas is known and 
the particle diffusion mode and energy losses are accurately 
known, whereas (ii) radio emission enables deduction of Up 
from electron synchrotron flux and spectrum once assump- 
tions have been made on the link between energetic electrons 
and protons, and between these particles and the magnetic 
field; and (Hi) assuming a SN origin for CR, SN statistics 
for a given (region of a) galaxy leads to an estimate of Up 
there. 

A substantial agreement among estimates based on the 
three methods is reached for most of the galaxies in Ta- 
ble 1. The only exceptions are (a) the SMC and (b) the 
SB/Seyfert-2 galaxy NGC 1068. Possible reasons for these 
discrepancies are: 

(a) The SMC mean magnetic field might be sufficiently weak 
that there is only a partial CR confinement, so that most CR 
diffuse out to the intergalactic space (Abdo et al. 2010d). If 
so, the radio method yields the (lower) actual particle energy 
density, whereas the SN method gives an estimate of the 
(higher) produced amount. 

(b) Unlike for NGC 4945 (the other SB/Seyfert-2 galaxy in 
our sample), the spectral energy distribution of NGC 1068 
differs markedly from that of purely SB galaxies, suggesting 
that its 7-ray emission may be unrelated to star formation 
(Lenain et al. 2010). Indications that this may be the case 
are that (i) there is no galactic wind emanating from the SB 
region of NGC 1068 (unlike all other SB galaxies in Table 1), 
and (ii) its location in the diagnostic /12//25 vs. /60//100 
plane 1 (see Persic & Rephaeli 2007 and references therein) 
suggests that NGC 1068 does not behave as a typical SB 
galaxy. (Its colors suggest an 'evolved starburst phase', in 
which star formation activity has decreased dramatically 
from the SB peak, and cirrus emission, which is unrelated 
to ongoing star formation activity, apparently contributes 
strongly to the IR emission.) Finally, NGC 1068 does not 
conform with the trend set by the other star-forming galax- 
ies in Table 1, in the L(>100MeV) vs. ^sNxM gas plane 
(Fig. 3d) and the L(>100MeV) vs. i/ S n plane (Fig. 3e; see 
also Lenain et al. 2010). The discrepancy is even larger if 
the true ^sn is significantly lower then the value reported 
in Table 1 and used in Fig. 3 (as it would be if derived from 
the SB-only component of the circumnuclear IR emission). 

Once these considerations are weighed in, the discrep- 
ancy between the very high CRp energy densities (U p ~ 
O(10 2 ) eVcm -3 ), deduced for SB nuclei, and the low values 
(U p ~ O(10 _1 )eV cm ' ), deduced for very quiet environ- 
ments, appears very significant. 

Observationally, the link between HE 7-ray emission 
and SNe is displayed in Fig. 3 for our sample galaxies. The 



1 The IRAS flux densities of NGC 1068 at 12, 25, 60 and lOO^m 
arc /12 = 38.7, / 28 = 87.4, f 60 = 190, /100 = 277. 
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Figure 3. Correlations for all the star-forming galaxies detected by Fermi/LAT and reported in Table 1 (empty circles: star burst /Seyfert- 
2 galaxies; filled squares: starburst galaxies; filled triangles: quiescent galaxies. Arrows indicate HE7-ray flux upper limits: Arp220 and 
M33 at, respectively, higher and lower flux limits.) The Milky Way luminosity has theoretical uncertainties stemming from the range 
of possible emission models matching the Galactic diffuse 7-ray background (sec Strong ct al. 2010). Wherever plotted, SN rates are 
displayed with a factor of 1.5 uncertainty (see text), except for NGC4945 and NGC 1068 where the actual range of deduced values is 
reported in Table 1. (a) Correlation between 7-ray luminosity, SN rate within the region of radius r s under study, CRp residency time, 
and total gas mass, (b) Distribution of gas mass as a function of luminosity, (c) Distribution of residency time as a function of source 
size. To guide the reader's eye, the dotted line drawn through the data represents the T rcs oc rf relationship, (d) Correlation between 
7-ray luminosity, SN rate, and total gas mass, (e) Correlation between 7-ray luminosity and SN rate. For comparison, the dotted line 
drawn through the data shows the uh^ slope, implied by the observational L oc SFR 14 relationship reported by Abdo et al. (2010b). 



theoretical link is apparent from combining Eqs.(9,13) 

L oc i^sn r ros r~ 3 M gas . (14) 

In deriving Eq.(14) we considered that the volume-averaged 
7-ray emissivities of the sample galaxies, appearing in 
Eq.(9), vary little from galaxy to galaxy. This is due to the 
relatively uniform GeV-TeV photon slopes, and correspond- 
ing CRp spectral slopes, observed for these galaxies (Abdo 
et al. 2011). Panel (a) shows the observational correlation 
between 7-ray luminosity, SN rate within r s , CRp residency 
time, and total gas mass; the plot shows the relation ex- 
pressed in Eq.(14). Can it be reduced to fewer parameters? 



It can, for example, if r ros oc rf; observationally, this latter 
correlation roughly holds for our sample (see panel (b)). A 
possible reason for this: if in most galaxies in the sample 
the CRp energy loss is dominated by p-p interactions, so 
Trcs ~ Tp P oc n _1 , and M gas is roughly constant within the 
sample, so n oc r^ 3 , then r ros oc r 3 follows. Indeed, in most 
of our galaxies the CRp energy losses are dominated by p-p 
interactions (see Table 2) and the gas mass varies relatively 
little (a decade) within the sample (see panel (c)). Panel (d) 
shows the 'reduced' relationship between 7-ray luminosity, 
SN rate, and total gas mass. A further step from Eq.(14), by 
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taking into account that M gaa is relatively constant across 
the sample, is represented by panel (e) that shows the 7-ray 
luminosity - SN rate relation; note that the distribution of 
data points is compatible with a i/g^ slope, which follows 
from vs~n oc SFR and the observational L 7 oc SFR 1 ' 4 rela- 
tion found by Abdo et al. (2010b) for this same sample of 
galaxies (except the composite SB/AGN sources). 

It should be noted that the steps leading from panel a 
through panel d to panel e are based on the relative con- 
stancy of M gas within the sample; however, this is a clear 
selection effect stemming from the fact that L(> 100 MeV) 
is higher for a larger M gas , because star formation is higher 
and because there is more target gas available for p-p inter- 
action. Given the low fluxes emitted by star-forming galax- 
ies, which are close to the sensitivities of the Fermi/LAT 
and Cherenkov instruments, it is clear that only the bright- 
est galaxies can be observed, i.e., the most gas-rich ones, 
with log Mm ~ 9 — 10 (see the HI mass function of galaxies, 
e.g. Zwaan et al. 2003). 

With the exception of (the mostly Seyfert-2 galaxy) 
NGC 1068, current measurements of the galaxies in our sam- 
ple provide appreciable evidence for the direct link between 
SN and VHE 7-ray emission. 



7 CONCLUSION 

Active star formation in galaxies leads to production of CR 
protons and electrons via the Fermi-I diffusive shock ac- 
celeration mechanism in SN remnants. Basic considerations 
suggest that the timescales required for CRp to be Fermi-I 
accelerated (r + ) and to lose energy via pion decay into pho- 
tons and e + e~ pairs or via advection by bulk outflows (r_), 
are shorter than the timescale of starburst activity in galax- 
ies (tsb); the latter is itself comparable with the dynamical 
timescale of galaxies. A consequence is that in a SB region a 
balance can roughly be achieved between energy gains and 
losses of CR during a typical burst of star formation. Un- 
der equilibrium conditions in a galaxy, a minimum-energy 
configuration of the field and the CR may be attained. This 
implies that energy densities of particles and magnetic fields 
may be in approximate equipartition. 

For 10 galaxies for which pointed GeV-TeV data are 
available (from Fermi/LAT) and from Cherenkov tele- 
scopes), we have estimated U p using the the radio, 7-ray, 
and SN methods. A substantial agreement among estimates 
based on the three methods is reached for 9 of the 11 galac- 
tic regions (3 SB nuclei, the Milky Way's disk and central 
region, 4 quiescent galaxies, 2 SB/Seyfert-2 galaxies) con- 
sidered in this paper. The discrepancy between the very 
high values (U p ~ O(10 2 ) eVcm" 3 ) found in SB nuclei 
such as M82, NGC 253 and NGC 4945, and the low values 
(U p ~ O(10 _1 )eV cm ) found in very quiet environments 
such as the Local Group galaxies, appears to be highly sig- 
nificant. 

The results of this study extend the validity of our pre- 
vious suggestions that: 

(i) CRp energy densities in star- forming galaxies range from 
O(10 2 ) eV cm~ 3 in very active environments (e.g., SB nuclei) 
down to O(10 _1 ) eVcm~ 3 in very quiet ones (e.g., Local 
Group galaxies); 

(ii) CR energy densities and magnetic fields, inferred from 



radio data under the assumption of energy equipartition, can 
be used as proxies of the actual quantities that are measured 
directly only from 7-rays. This could be particularly useful in 
the case of distant galaxies, whose (unbeamed) 7-ray fluxes 
are too faint to be measured; 

(Hi) core-collapse SN, both in quiescent galaxies and in SB 
galaxies, are likely to share a universal CR acceleration effi- 
ciency; and 

(iv) the Fermi-I acceleration mechanism, assumed to be at 
play in the environment of SN remnants, leads to accel- 
eration timescales for CRp in galaxies such that particles 
and fields attain equilibrium over typical star formation 
timescales, in agreement with observational evidence. 
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APPENDIX: PROTON/ELECTRON RATIOS 

Expanding on Schlickeiser (2002), let us assume that 
electrons and protons are accelerated out of a thermal 
bath, characterized by a non-relativistic kinetic energy 
To ~ lOkeV, to a differential PL spectrum (in momentum) 
Nj(p) — No,jP~ qi where j=e,p. At injection, q c — q p . Later, 
different loss and diffusion processes for the two CR species 
will in general result in q e / q p . If electrical neutrality of pri- 
mary CRs is conserved, then the total number of CR particle 
of each species per volume unit is 



f-oc poo 

i = / N c (T)dT = / N p (T)dT. 



(15) 



The energy-momentum relation T—y^m 2 c 4 + p 2 c 2 —mc 2 im- 
plies P=\/T 2 /c 2 + 2Tm and dp/dT=(T/c 2 + m)(T 2 /c 2 + 
2Tm)~ 1/2 . As Nj(T)=Nj[p(T)]dp/dT , the final expression 
is (see Schlickeiser 2002) 



N i( T ) = (T + mc 2 ) (^ + 2Tm) 



-(9j+l)/2 



(16) 



Inserting Eq.(16) in Eq.(15), we obtain the normalization of 
each CR species, 



No,i = n (-2- + 2T m) 



(9j-l)/2 



(17) 



Because of the assumed electrical neutrality of the primary 
CRs, we get 



N, 



o,p 



_1 [m p + (To/2c : 



2 VI (<?p-1)/2 



(2T 



,(<2p-l)/2 



iVo, c q c -l [ mc + (To/2 c 2 )]^- 1)/2 (2To)(--D/2 

gp-1 (2T m p )^- 1 ^ 2 
q e -l (2T m c )(9e-i)/2 

since To << m c c 2 << m p c 2 . 

• The p/e number density ratio is 

N P (T) dT 



(18) 



"(T; M ) = Ne{T)dT - 
Inserting Eqs.(16),(18) in Eq.(19), we obtain 



(19) 



oj(T;q p ,q e 



(? p -l) (2T m p c 2 ) — 



(*> -1) (2Tom e c 2 ) a£ 2^ 



T 1 -^ (T + m p c 2 )(r + 2m p c 2 )- 
T!-9e (T + m c c 2 )(T + 2m e c 2 )~ 



9p + l 



(20) 



From Eq.(20) it is straighforward to see that at injection, 
when q p = q c = g in j, it is (Schlickeiser 2002) 



u(T;q ini ) { 



= 1.. 

oc 



\ m p c z / 

— (^p) 2 
V rn c J 



T/c 2 «m c 

m c « T/c 2 «m p (21) 



T » m p c 2 



Notice that for <? inj ~ 2.2, cj(T >> 1 GeV) ~ 10 2 , as re- 
marked by Schlickeiser (2002). 

• The p/e energy density ratio is 

J™N p (T)TdT 



n(q P ,q a ) 



CN e (T)TdT 



(22) 
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Inserting Eqs.(16),(18) in Eq.(22), we obtain 



(«e-l) (2T m e c 2 ) 2 ^ 

/ T 2 (T+m p c )(T+2m p c ) 2 dT 

x "".r 1 o o — ' *- 23 - ) 

T T _ (T+m c c 2 )(T+2m e c ! ) "T - dT 

JT 

Calling the integrand functions on the top and bottom of 
the r.h.s. of Eq.(22), respectively, f P (T) and f e (T), we can 

rewrite /~ / p (T)dT = / p (T)dT+/~ c3 f p {T)AT and 

/~ / e (T)dT = J™^ / e (T)dT + JJ oc2 / e (T)dT. Hence, in 
the approximations T « m p c 2 ,T >> m p c 2 and, respec- 
tively, T « m c c 2 , T >> m e c 2 , we can write Eq.(23) as 



(2Tom e c 2 ) 2 ^ 

Qp — 1 iTTlpC 2 IOO 

(2m p c 2 ) 2— 1^T 3 -«P + ^J^T 2 -"P\ 

x — — ^ It " , 2 ' 9P , (24) 

9c - 1 , |TT1 C C^ , loo ' V / 

(2m c c 2 ) — 3— l_T3-9o + l_T2-?o 

where T = 10~ 5 GeV, m c c 2 = 0.511 x 1(T 3 GeV, m p c 2 = 
0.94 GeV. Note that at injection q p — g c = qi n j, so that 

K (g inj ) ~ (!^E)^- ... T»lGeV. (25) 

Eq.(25) is the counterpart of the third of Eqs.(21). 

Eq.(24) enforces electric neutrality of the primary CRs 
for any pair of q p , q B . In computing k from Eq.(24) we as- 
sume ye to be measured from radio data according to stan- 
dard synchrotron theory, q c = 2a + 1. However, the CRp 
spectral index, q p , is not directly measurable hence some 
guess must be made: 

(a) in starburst regions, that are sources of primary CRs, 
we assume the average CRp spectral index to stay close to 
the injection index, q p ~ qi n j-. as VHE7-ray observations of 
SN remnants suggest <j>i n j ~ 2.2 (Aharonian et al. 2006), we 
assume q p ~ 2.2 in starburst cores; 

(b) in quiescent galaxies like the Milky Way, small source re- 
gions are mixed with large passive regions, so that in the re- 
gions typically sampled in these galaxies, CRp are mostly af- 
fected by confinement effects only (leaky-box model: homo- 
geneous distribution of CRp sources, no CRp energy losses, 
CRp confinement time scaling inversely with energy), and 
the steady-state solution of the diffusion-loss equation yields 
<7 P = <7inj + S 2 . Because S ~ 0.5 (Gupta & Webber 1989), 
we assume a spectral index q p = 2.7 for CRp in quiescent 
galaxies. 



2 The general diffusion-loss equation is dN g^ = DV 2 N(E) — 
-^fi + JL[b(E)N{E)\ + Q(E) (Ginzburg & Syrovatskii 1964) 
where D is the scalar diffusion coefficient, Q(E) is the source 
term, r(E) is the energy-dependent time confinement timcscalc, 
b(E) is the energy loss rate, N(E) is the steady-state distri- 
bution of particles. At steady state (^rzrjfO, for a homoge- 
neous distribution of sources (DV 2 N(E) = 0) and negligible 
CR energy losses (b(E) = 0), the loss-diffusion equation becomes 
N(E) = Q{E)t(E) (leaky-box model). Since t(E) oc E~ s , as- 
suming Q(E) oc £-9i„j we get AT(_E) ^ £-(<7i„j+<5). 



